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ABSTRACT 

We present a study of ten i?-meson decays to a D^*\ a proton-antiproton pair, and a system of 
up to two pions using BABAR's data set of 455 x 10^ BB pairs. Four of the modes {B — > D^pp, 
B — > D*^pp, B — > D^ppK~ , B — > D"'^^'^pp-K~) are studied with improved statistics compared 
to previous measurements; six of the modes {B~ D^ppn^ , B^ ^ D*^pp'K~ , B ^ D^pp7r~7r~^ , 
B D*^pp'K~7r^ , B~ ^ D'^ppir^TT^ , 5^ — > L)"**^pp7r~7r~) are first observations. The branching 
fractions for 3- and 5-body decays are suppressed compared to 4-body decays. Kinematic 
distributions for 3-body decays show non-overlapping threshold enhancements in m{pp) and 
m{D^'^p) in the Dalitz plots. For 4-body decays, m{pTT~) mass projections show a narrow peak 
with mass and fuh width of (1497.4 ± 3.0 ± 0.9) MeV/c^ and (47 ± 12 ± 4) MeV/c^, respectively, 
where the first (second) errors are statistical (systematic). For 5-body decays, mass projections 
are similar to phase space expectations. All results are preliminary. 
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I. INTRODUCTION 



The decays' of B mesons to final states with baryons have been explored much less systematically 
than decays to meson-only final states. Such decays have their own distinctive features; in par- 
ticular, the suppression of the rate for two-body decays and the rate enhancement for low masses 
of the baryon-antibaryon system in multi-body decays [1-7]. The first observed exclusive decays 
were the CLEO measurements of B ^ A'^ pTr{n) [8] followed by B^^D*^ppTr~ and B^ ^ D*~pn^ 
[9] supporting a prediction [10] that the decays with Ac are not the only significant contributions 
to the baryonic B decay rate and that A^A^'+anything is also important, where A^ and A^' are 
nucleons. Belle observed B^ ^D^^pp [11] obtaining a branching fraction about five times smaller 
than for D*~^ppTr~ . 

We report the branching fractions and kinematic distributions for ten baryonic B decays: 



3- body decays B^ - 

4- body decays B^ - 

B- 

5- body decays B^ - 

" B- 



• D^pp 

•D^ppTT^ 
^D^ppTT^ 

• D^ppTT~TT^ 
D'^ppTT^TT'^ 



and D*^pp, 
and D*'^ppTT~, 
and D*^pp-K~ , 
and D*^ppTT~7r~^ 
and D*~^ppTT~'ir' 



(1) 



where the latter six in the list are first observations.'^ The D* are reconstructed as D^tt~^ and 
D^TT^; D as K~7r~^, K~tt~^tt^, K~ tt'^ 'k~ tt'^ , and K~tt~^it'^ . This makes up 26 reconstructed decay 
chains (3 for each of eight B decays with and 1 for two with D~^). This study uses 455 x 10^ 
BB pairs and supersedes the previous BABAR publication of B^ ^D^^^pp and D^^ppir^ [13] using 
232 X 10^ BB pairs. Figure 1 shows the typical valence-quark diagrams for 3- and 4-body decays. 

Until recently, interest has focused on the dynamical features of baryonic decays, on studying 
rare modes (both 6 — > u and 6 — > s transitions), and on the role these modes play in accounting 
for the overall production of charm in B decays. Many theoretical studies have appeared [10, 14- 
27] using the decays as a playground for phenomenological models including some with pp bound 
states [25] and multi-quark intermediate resonances [26]. We hope to shed light on these models 




^(*)o 



B- 




(a) Diagram for B°^D'-*'>"pp 



(b) Diagram for B —*D^*^"ppTv 



FIG. 1: Typical valence-quark diagrams for (a) 3- and (b) 4-body B decays. 



' Charge conjugation of particles and decays are implied throughout this document unless otherwise stated. 

The charge conjug ate B^^D* pn is not included since CLEO only detects antineutrons n. 
^ B^ D^*^''pp are suppressed by [12] with respect to Fig. la for b^uD'^*''^ and are beyond our sensitivity. 
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by studying the relatively unexplored territory of baryonic B decays involving a D^*^ meson. 

II. THE BaBar detector AND DATA SET 

The BABAR detector is described in detail elsewhere [28]. Exclusive i?-meson decays are simulated 
with the Monte Carlo (MC) event generator EvtGen [29] and hadronization is simulated with 
Jetset 7.4 [30]. We use GEANT4 [31] to model interactions of particles traversing the detector, 
taking into account the varying detector conditions and beam backgrounds. 

This study uses the data set of 414 fb~^ collected with the BABAR detector at the PEP-II 
asymmetric-energy e~^e~ collider at the SLAC National Accelerator Laboratory. Figure 2 is a data 
event display of B^ ^ D^pp followed by ^ K~^7r~ . For each daughter particle, we can see the 
drift chamber hits matching up with projected Cherenkov light cones in the ring-imaging detector. 

III. ANALYSIS METHOD 
A. Event selection 

We select events with a B candidate reconstructed in one of 26 decay chains in Eqn. 1 in two steps: 
First, the BABAR data set is reduced by requiring the presence of a p, p, and D. The average 
momentum of a proton produced in a typical B decay listed in Eqn. 1 is 1 GeV/ c. Protons are 
identified with a likelihood-based selector using information from the silicon vertex tracker, drift 




FIG. 2: Event display of B^ ^ D^pp, ^ K'^ir transverse to the e beam direction. 



9 



chamber, and ring-imaging Cherenkov detector, which has a 98% selection efficiency and a 1% kaon 
fake rate. The D mesons decay to charged tracks and vr^s. Charged tracks are required to be in the 
fiducial volume and have a distance of closest approach to the beam spot < 1.5 cm. Neutral pions 
are formed from two well-separated photons with 115 <m^^ < 150MeV/c^ or from two unseparated 
photons using the second moment of the electromagnetic calorimeter energy distribution. The 
same requirements apply to the non-composite daughters of B and D* decays. 

Second, the data set is further reduced by the explicit reconstruction of B candidates. B 
candidates are formed with a D mass [32] constraint using a Kalman fitter [33] with its vertex 

probability >0.1%. To suppress continuum e+e~ ^ qq {q = u,d, s,c) events, we compute the 
angle between the thrust axes of the B candidate and the rest of the event [34] . For modes with a 
D*, the mass difference Am = Tn(D^TT) — m{D^) is required to be within 3a of the nominal value, 
where the resolution is around 0.8MeV/c^. For B^ ^ D'^ppir" it~^ , we require Am> 160MeV/c^ on 
the D^7:~^ system. For all decay modes, D candidates are required to be within 3a around the 
nominal value, where the mass resolution for Kn, Kmr^, KninT, and Kmr is around 6, 10, 5, and 
5 MeV/c^, respectively. Furthermore, the average momentum of a kaon produced in a typical D 
decay is 0.9GeV/c and is found with the above-mentioned likehood technique, which has an 85% 
selection efficiency and a 2% pion fake rate. For K-ktt^ decays, we use the squared decay amplitude 
based on a Dalitz plot model [35]. Lastly, requirements are optimized by maximizing the ratio of 
the squared expected signal yield and its sum with expected backgrounds. For cut values with 
broad maxima, they are chosen to be uniform across related decay modes. 

After events are filtered through the two steps, we are left with 0(10^) for all modes. The average 
number of B candidates per event that pass all requirements ranges from 1.0-1.7, increasing with 
the multiplicity of the decay. If there are multiple such candidates, then we choose one with 
the L>(*) mass closest to the nominal value; furthermore, if candidates share a Z)", we choose 
one at random. The reconstruction efficiencies are found using MC containing the desired B 
decay. In general, efficiency decreases with particle multiplicity: D^pp, Kir is highest at 19% and 
D*ppTnT, Kinnr is lowest at 0(1%); they are given later in Table Ila. 



where ^/s is the e^e~ center-of-mass (cms) energy and (Eb) is the B candidate momentum 
(energy) in the cms. Figure 3 gives ttt-es-^-^' scatter plots for the six B decays that are first 
observations. A concentration of correctly reconstructed B candidates is visible in the region 
of rriEs and AE near the nominal B mass of 5.28GeV/c^ and zero, respectively. The uniform 
distribution of dots over the entire plane is indicative of the general smoothness of the background 
events. 

A 2-dimensional (2d) probability distribution function (pdf) is used to fit m-Es-AE via the 
unbinned extended maximum likelihood technique [36]. The 2d pdf is a sum of two components 
for all B decay modes except for D*^pp'K, which requires three. These are denoted as P5, Pb, and 
Pp respectively, for the signal, background, and peaking background events. The motivation and 
description for Pp is given later and continued in Sec. IIIC. The parameters (yields) associated 
with the components are O5, 0^, and Qp {ns, ns, and np), respectively. Since the fit variables 
are uncorrelated to a good approximation, each 2d pdf component is written as product of two 
Id pdfs. Ps is the product of two functions with a Gaussian core and a power-law tail [37] 
written as Pg{m-Es',^s)Pg{AE;^}'g). Pp is the product of a threshold function vanishing at the 



B. Fit method and yields 



The signal B yield is extracted by fitting the joint distribution of 




and 




(2) 
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nominal B mass for mss and a 2 -order Chebyshev polynomial of the first kind for AE written 
as Pt{mEs;^B)PciAE]n'j^). Pp for D*^ppiT is written as Pg{mEs;^p)Pg{^E;n'p). 

The fit window of 5.22 < mss < 5.30GeV/c^ and \AE\ < 50MeV provides ample sideband 
regions because the signal resolution of rngs and AE are relatively narrow at 2.2-2.5 MeV/c^ and 
8-10 MeV, respectively. Signal slice projection plots for tties are defined to be within 2.5cr of the 
mean value of AE and sideband to be outside 4cr and vice versa for AE plots. As an illustration 
of the fit projection over the entire window, Fig. 4ab(cd) gives projections of the 2d pdf in the 
signal (sideband) slices. For ^ D'ppTnr, corresponding to the scatter plot in Fig. 3c, the good 
description of background in the sideband slices gives us confidence that it is well modeled in the 
signal box. 

The likelihood for N data events is defined as 

-(ns+ns) ^ <=^™t'^ 

C = n Pgi^ESi-, ^s) PgiAEf, n's) +nB PtimESi] ^b) PdAEf, Q'b) ] , (3) 

i = l ' V ' ^ V ' 

signal pdf Pg background pdf Pg 
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FIG. 3: m-ES-^E scatter plots for new observations of B decays: (a) B 
B^^D^pp'K-n, (c) B-^D+pp-jnr, (d) B-^D*^pp'K, (e) B^ ^ D*^ ppin: , and (f) B' ^ D*+ pp'K-K . 
We show the cleanest mode with D^'~^ ^K~tt~ , K~Tr~^7r~^, resp. Histograms on top (right) are 
mES i^E) projections in 1 MeV/c^ ( MeV) bins; no cut is made on the complementary variable. 
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where ruESi ™d AEi are the values for the i event. For D ppn, we add a thhd component in 
Eqn. 3 for peaking background. The quantity — ln£ is minimized with respect to the parameters 
and yields using Minuit [38] via the RooFiT toolkit [39] in ROOT [40]. We fix some parameters 
to values obtained in fits to MC distributions. For Ps, these are the Gaussian width and power-law 
tail parameters; for Pb, the end-point of the threshold pdf for rfiEs! and for Pp, all parameters 
including the the yield np using the branching fraction and acceptance found for D*^pp7r. To 
check for possible pathological behaviors in the fit, we perform one thousand mock experiments 
where events are drawn from the pdf and subsequently fit with the same pdf. All means (widths) 
of signal yield pull distributions are consistent with zero (unity). 

Figs. 5, 6, and 7 give the collection of tties sliced projections in the AE signal region for 
3-, 4-, and 5-body B decays, respectively. Signal yields range from 50-3500 events for the ten 
B decays; they are given later in Table Ila. Lastly, additional work is done for four B^ decay 
chains D^*' ppinT, KniT^ and £>(*)° ppiTTT, KirTTTT shown in Fig. 7bcef because the fits converged to 
non-physical Gaussian mean values for Pg{m-Es) and Pg{AE) due to the low signal-to-background 
ratio. For these, the means are fixed to the values found for the Kir counterpart and the end- 
point for Pt{mEs) is floated. We note that these measurements do not contribute significantly to 
the branching fraction averaged over decay modes because of large systematic uncertainties 
compared to the corresponding Ktt measurement; this will be discussed later in Sec. IIIC. 

C. Systematic studies 

Table la lists the individual sources of systematic uncertainties for the branching fraction mea- 
surement in five categories (I-V): B counting, assumed branching fractions for T(4S') and D^*\ 
acceptance, pdf, and peaking background estimate. 

I. BB is counted by subtracting the total number of hadronic events taken at ^/s of 10.58 GeV by 
the expected number due to continuum events, which is estimated using data taken at 10.54 GeV. 
The total uncertainty is 1.1% and the largest component comes from the comparison of detection 
efficiencies of hadronic events in data and MC. 

II. We assume equal production of T(45) — > B^B^ and B^B~ as well as nominal branching 
fractions for D* and D. (1) For T(45), we take the absolute difference between the nominal value 
and 50%, which gives 1.6%. (2) For D*^ (£>*+), the uncertainties are 4.7% (0.7%). For D, the 
uncertainties are is 1.3%, 3.7%, 2.5%, and 2.3% for Ktt, Ktt-k^, Kittt-k, and Ktttt, respectively. 




mEs(GeV/c^) AE(GeV) mEs(GeV/c^) AE (GeV) 



(a) -18.8<A£;<14.7 (b) 5273.6 < mES < 5284.8 (c) A£;<-28.9 or > 24.8 (d) mES<-5270.3 

FIG. 4: The m-Es-AE fit projections for B~ ^ D^ppinT, D'^ ^ Ktttt in (ab) signal and (cd) 
sideband regions with requirements denoted in the caption in MeV/c^ (MeV) for tties i^^)- 
signal slices in (ab), the top blue curve is the signal pdf Ps] bottom red is the background pdf 
Pb- For sideband slices in (cd), no part of the signal pdf is present in this region and the 
background pdf is identical to the corresponding plot in (ab). 
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III. The reconstruction efficiency has seven contributing sources. (1) Charged track identifica- 
tion uncertainty of 0.5% per track is obtained by comparing data and simulations with e~^e~ ^t~^t~ 
events where one tau decays leptonically and the other tau decays hadronically. (2) Because soft 
charged pions from 1)*+ decays leave hits only in the silicon vertex tracker, they have an addi- 
tional uncertainty of 3.1%, which is obtained from using the helicity angle distribution of the decay 
daughters [41]. (3) Neutral pions identification of 3.0% is obtained by comparing data and MC 
r+T~ events where one tau decays leptonically and the other decay includes a vr*^. (4) Acceptance 
of signal B decays computed using the MC assuming the uniform phase space decay model is com- 
pared with the data distribution in bins of m?{pp) and rin?{D^p) where the most variation is seen. 
We obtain uncertainties ranging from 0.8-9.7% with modes containing a D*^ at the higher end of 
the spectrum. (5) Kaon identification uncertainty of 0.5% is obtained by comparing data control 
samples of D^^Ktt from inclusively produced D*^D^7r and from similar B decays. (6) Proton 
identification uncertainty of 1% is obtained by comparing various data control samples of inclu- 
sively produced A^pvr" with decay topologies of high charged-particle multiplicities that mimic 
our B decay environment. (7) Corrections made for particle identification contribute additional 
statistical uncertainties that range from 1.5-2.5%. 

IV. The fit pdf has four contributing sources. (1) As described in Sec. Ill, some parameters 
are fixed to values determined by the MC sample. These are varied by the uncertainties obtained 
when fitting the MC and the quadrature sum of the typical yield changes are 1.3, 2.8, 5.7, and 
3.4%for modes with D — > Kit, Ktttt^, KmriT, and Ktttt, respectively. (2) A more general signal 
pdf is chosen and the typical yield change is 0.6%. (3) A more general background pdf is chosen 
and the typical yield changes are 0.8, 4.5, 1.3, and 2.0%for modes with D^Kn, Kmr^, K-k-rtt, 
and Ktttt, respectively. (4) Sec. IIIB described a thousand mock experiments to check for possible 
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pathologies of the pdf. We generahze this study by embedding the MC with signal B decays into 
background events drawn from the background pdf. Signal yield bias ranges from 0.4-2.2%, which 
we assign as a systematic uncertainty. 

V. Peaking background sources have three contributing sources. (1) Possible sources of events 
that peak in the signal box around uies = 5.28GeV/c^ and AE = OM.eV have been studied exten- 
sively using a sample of O(IO^) MC events with B decays hadronized using Jetset 7.4 [30] as well 
as a sample of 0(10®) MC events dedicated to a set of B decays to a D^*\ a pp pair, and a system 
of up to two charged or neutral pions followed by a comprehensive set of D^*) decays including 
those not reconstructed in our study. We find only one scenario meeting our requirement: events 
generated as B^ —>■ D*~^pp7:~ reconstructed as B~ D*^ppiT~ . In this case, we fit the MC and 
add the pdf component in the fit to data in Eqn. 3 as described in Sec. IIIB. The normalization 

TABLE I: Systematic uncertainties for measurements of B branching fractions, B. 
(a) List of individual systematic uncertainties in five categories (I-V). 



Item Source description %of B 



I 


BB counting 


1.1 


II. 1 


Branching fraction of T(4S') 


1.6 


II.2 


Branching fraction of D^Ktt, Ktitt^ , Ktttttt, Ktttt 


1.8, 4.4, 3.2, 3.6, rcsp. 


II.3 


Branching fraction of D* — ^D^tt^, D°7r+ 


4.7, 0.7, resp. 


III.l 


Efficiency of finding charged tracks not including soft pions 


0.5 per track 


III.2 


Efficiency of finding soft charged pions from 


3.1 per TT 


III.3 


Efficiency of finding neutral pions 


3.0 per 7r° 


III.4 


Efficiency of finding B decays in bins of m?[pp) and m?[D^*''p) 


0.8 to 9.7 


III.5 


Efficiency of finding kaons in B decays 


0.5 per K 


III.6 


Efficiency of finding protons in B decays 


1.0 per p 


III.7 


Efficiency of particle identifcation based on data control samples 


1.5 to 2.5 


IV. 1 


Pdf parameter variation for modes with D^Ktt, Ktttt'^, Ktttttt, Ktttt 


1.3, 2.8, 5.7, 3.4, resp. 


IV. 2 


Pdf choice for signal events 


0.6 


IV.3 


Pdf choice for background events for Z)— s-A'tt, Kittt'^, Ktttttt, Ktttt 


0.8, 4.5, 1.3, 2.0, rcsp. 


IV.4 


Pdf yield bias by fitting mock experiments embedded with MC 


0.4 to 2.2 


V.l 


Peaking background in m-Es-AE 


8.0 for D*°ppTT 


V.2 


Peaking background in mES only {D^*^^pp7TTr, Ktttt^ and D^*>^ppTnr, Ktttttt) 


0.0 to 14.5 (77 to 85) 


V.3 


Peaking background from identical final states without a D meson 


0.5 to 13.5 



(b) Uncorrelated sources (%of B) (c) Correlated sources excluding II. 2, II. 3 (%of B) 



iV-body 


D decay 

B decay 


Ktt Ktttt'^ K 


TTTTTT 


Ktttt 


D decay 

B decay 


Ktt Ktttt'^ Ktttttt Ktttt 


3-body 




■^D°pp 


2.7 


5.5 


4.7 




B"- 


-^D°pp 


3.5 


6.9 


7.0 


II 




-*D*^pp 


2.2 


8.6 


8.7 






■^d*°pP 


4.6 


8.6 


7.7 


4-body 




^ D^ppTT^ 








5.7 


S"- 


^D^ppTT~ 






- 5.6 


II 


W- 


D*^ ppTT^ 


4.2 


7.6 


9.9 




5°- 


^ D*^ ppTT~ 


6.3 


8.6 


8.9 


II 


B-' 


D^ppTT^ 


14.5 


81.2 


77.3 




B-- 


D'^ppTT^ 


3.9 


7.0 


7.0 


II 


B-- 


D*^ppTT^ 


13.8 


86.3 


85.4 




B-- 


-^D*"ppTT- 


4.9 


8.8 


7.7 


5-body 




D'^ppTT^TT^ 


4.4 


8.8 


11.6 




B^- 


-^D^ppTT^TT^ 


4.2 


7.2 


7.3 


II 




^ D*'^ppTT^TT^ 


10.4 


11.6 


13.2 




B°- 


-^D*^ppTT^TT^ 


5.2 


8.9 


7.9 


II 


B-- 


->■ D^ppTT^TT^ 








15.0 


B-- 


ppTT^ TT^ 






- 5.9 


II 


B-- 


D*^ PpTT~ TT~ 


5.9 


14.9 


19.0 




B-- 


D*'^ PpTT~ TT~ 


6.8 


9.0 


9.2 



16 



is the product of the acceptance and branching fraction for D*~^ppTT~ . The total uncertainty on 
the latter quantity is 8% and is assigned as a systematic. (2) Possible sources of events that peak 
either in ttt-es = 5.28 GeV/c^ or in A£' = OMeV have been studied using the above-mentioned set 
of MC. We find that only a few samples exhibit a peak in tties that occur when a related decay 
is misreconstructed by losing or gaining either a vr'' or a 7 from the other B decay in the event; 
we find no cases of events peaking only in AE. In the rriES-peaking cases, the distributions peak 
broadly with a Gaussian width of around 10 MeV/ and has a smooth variation in AE. One such 
distribution is fit and the component is added to our nominal fit in Eqn. 3 as a systematic study. 
The signal yield changes varied from 0-15% with four exceptions: in the decays D^*^ppinr, Ktttt^ 
and D^*^ppTT7T, Rtttttt shown in Fig. 6bc (6ef), the signal yields varied by 77-85%. This is not 
surprising because we saw in Sec. IIIB that the fit was not stable against our initial strategy and 
had to be forced to find a peak at the desired location by fixing the means of Ps- For these four 
B modes, the Ktt measurement is heavily weighted due to these systematic errors with respect 
to both Knir'^ and KmnT when we average the B branching fraction over D decay modes. (3) 
Possible sources of events that share the final states with desired decays are studied by using m{D) 
sidebands in data. One such example scenario would be B ^ KcP'K^ , Kc^pKir misreconstructed 
as B^^D^pp, ^ Ktttt^ , which would live under the D mass peak in the m{D) distribution. To 
correct for such a bias, we scale the B signal yield in the m{D) sideband region to the amount of 
D background in the m{D) signal region. This is given later as a correction factor Uf, in Table Ila 
and the uncertainties range from 0.5-13.5%. 

Combining systematic uncertainties requires us to divide uncorrelated sources from those that 
are not. Table lb gives the uncorrelated items III. 4, III. 7, IV. 4, V.l, V.2, and V.3 summed in 
quadrature. Table Ic gives the correlated sources, which are computed in two steps. Correlated 
sources for the individual decay chains — III.l, 2, 3, and 6 — are summed linearly; Correlated sources 
for B decays — I; II. 1; III. 5 and 6; and IV. 1, 2, and 3 — are summed in quadrature. Lastly, the error 
matrix, V, is constructed. For decays with a D^, V is a 3x3 matrix spanned by the three decay 
modes {Ktt, Ktttt^, KirmT) labeled by a and /3; for those with a D~^, it is a 1 x 1 matrix for Ktttt: 

Vsyst 
. " ^ 

V = Vstat + V 

unc 

+ V 

cor is the sum of statistical and systematic components 

Vstat = diag[(T|^] for statistical errors from fitting for the signal yield (4) 

Vunc = cliag[(7^nj. ^] for uncorrelated systematic errors 

Vcor = diag[(T^oj, ^] + Pa/3 ^cor.a Ccor,/3 for Correlated systematic errors, 

where pap is the correlation coefficient and as is given later in Table Ila. The pa/s for branching 
fraction are given in [32] and all other coefficients are assumed to be unity. 



IV. BRANCHING FRACTIONS 
A. Weighted average 

Branching fractions are found in two steps: they are computed for 26 decay chains, then averaged 
over D decay modes. The branching fraction for a B decay followed by D^Ktt denoted by a, is 

where N = x 10^, B{D^*^) is the D^*) branching fraction, e is the acceptance, ns is the signal 
yield, and is the D-mass-sideband peaking background correction. Branching fractions for eight 
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B decays with a form a vector B = {Br-k, BxttttO^Bk-ktv-k)'^] for two with a B = {Bk-ktv)'^- 

Table II presents the branching fractions and their ingredients for the 26 decay chains and the 
combined values, B^^^g = B • w, which are averaged over decay modes using a method [42] to 

TABLE II: B branching fractions, B, (a) for 26 decay chains and (b) averaged over D modes, 
(a) B and its ingredients. See Eqn. 5. 



A^-body B decay chain 


ns±as 




e (%) B 




D' (%) 




3-body 




^D^pp, Ktt 


351±20 


7.6 


19.04 


3.89 




1.02±0.06 


II 


B^- 


-> D°pp, Ktttt^ 


431± P 


23.7 


7.01 


13.5 


- 


0.95±0.06 


II 


B^- 


-y D^pp, Kmnr 


448± 


10.1 


9.85 


8.10 




1.21±0.07 


II 


B°- 


-^D*°pp, D^n, Kn 


110± 


-1.4 


9.40 


3.89 


61.9 


1.08±0.12 


II 


B°- 


^D*"pp, D^TT, Ktttt" 


148±15 


3.9 


3.24 


13.5 


61.9 


1.17±0.12 


II 


B^- 


^D*°pp, D^TT, Ktttttt 


95± \i 


5.5 


5.15 


8.10 


61.9 


0.76±0.12 


4-body 


B°- 


■^D'^ppTT, Kttt: 


1816± 11 


55.2 


12.64 


9.22 


- 


3.32±0.10 


II 


B°- 


■^D*+ppTT, D°7T, Ktt 


392zb 20 


2.3 


6.79 


3.89 


67.7 


4.79±0.26 


II 


B°- 


■^D*+ppTT, D^TT, KtTTT^ 


601± 


20.7 


3.08 


13.5 


67.7 


4.53±0.22 


II 


B"- 


^D*^ppTr, D'^tt, Kmnr 


378± 11 


19.9 


3.66 


8.10 


67.7 


3.92±0.24 


II 


B-- 


-^D'^ppTT, Kit 


1078± 


13.1 


15.89 


3.89 


- 


3.79±0.14 


II 


B-- 


-^D^ppTT, Knir^ 


1176± II 


41.1 


5.53 


13.5 




3.34±0.16 


II 


B-- 


-^D^ppTT, Ktttttt 


1296± % 


33.0 


7.82 


8.10 




4.38±0.20 


II 


B-- 


-^D*°ppTT, D°Tr, Ktt 


328±22 


2.1 


7.71 


3.89 


61.9 


3.86±0.26 


II 


B-- 


-^D*°ppTT, D°TT, Ktttt^ 


482± 11 


46.5 


2.87 


13.5 


61.9 


3.99±0.32 


II 


B-- 


D*'^ppTT, D^TT, Ktttttt 


343± 11 


32.4 


4.04 


8.10 


61.9 


3.37±0.34 


5-body 


B°^ 


^ D^ppTTTT, Ktt 


438± 11 


7.7 


8.19 


3.89 




2.97±0.22 


II 


B°- 


^D'^ppTTTT, Ktttt^ 


663± 11 


155.2 


2.92 


13.5 




2.83±0.36 


II 




^ D^ppTTTT , Ktttttt 


770± 


39.7 


3.75 


8.10 




5.28±0.48 


II 


B°^ 


^D*"ppTnT, D°TT, Ktt 


61±12 


1.8 


2.89 


3.89 


61.9 


1.87±0.38 


II 


B°^ 


^D*°ppTnT, D°TT, Ktttt" 


142± 33 


36.7 


1.27 


13.5 


61.9 


2.19±0.66 


II 


B°- 


^D*^ppTTTT, D^TT, Ktttttt 


163± 11 


12.8 


1.33 


8.10 


61.9 


4.93±0.99 


II 


B-- 


-^D~^ppTTTT, Ktttt 


475± 11 


6.6 


6.74 


9.22 




1.66±0.13 


II 


B-- 


D*^ppTTTT, D^TT, Ktt 


57± 9 


-12.4 


2.93 


3.89 


67.7 


1.98±0.26 


II 


B-- 


D*^ppTTTT, D^TT, KtTTT^ 


94± \i 


-0.6 


1.25 


13.5 


67.7 


1.82±0.27 


II 


B-- 


D*^ppTTTT, D^TT, Ktttttt 


66± 


4.8 


1.52 


8.10 


67.7 


1.61±0.32 



(10-4) 



(b) B averaged over D decay modes and previous measurements [9, 11, 13] in units of (10 See Eqn. 6. 



iV-body B decay B±aste.t±asyst P{x^) Rcfs. [9, 11] Ref. [13] 



3-body 


B°- 


^D^pp 


1.02±0.04±0.05 


4.3 


12% 


1.18±0.15±0.16 [11] 1.13±0.06±0.08 


II 


B"- 


^d*"pP 


0.97±0.07±0.09 


4.1 


13% 


1.20±g:i ±0.21 [11] 1.01±0.10±0.09 


4-body 


B"- 


^ D^ppTT^ 


3.32±0.10±0.27 






3.38±0.14±0.29 


II 


B°- 


D*^ ppTT^ 


4.55±0.16±0.37 


1.2 


54% 


6.5 ±\l ±1.0 [9] 4.81±0.22±0.44 


II 


B-- 


D^PpTT^ 


3.72±0.11±0.23 


3.4 


19% 




II 


B- 


D*^ppTT^ 


3.73±0.17±0.39 


0.5 


79% 




5-body 


B°- 


-> D^ppTT^TT^ 


2.99±0.21±0.44 


0.3 


85% 




II 


B"- 


D*'^ppTT^TT^ 


1.91±0.36±0.29 


0.5 


78% 




II 


B- 


D^ppTT^ TT^ 


1.66±0.13±0.27 








II 


B- 


D*'^ PpTT~ TT~ 


1.86±0.16±0.18 


0.2 


91% 
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compute the weights from error matrices in Eqn. 4 



V 



u 



w 



(6) 



where u = (1,1,1)"^. The squared combined statistical uncertainty is cTg^at = w'^-^sta.fW and similarly 

computed for systematics. All measurements are statistically significant. 

For B decays with a meson, consistency between the three measurements is evaluated by 
= iB — Bnvgu)'^ ■ • {B — Bs.vgu) and the associated probability P(x^) = exp[— x^/2]/2. For 

all cases, P(x^) is greater than 10%. 

B. Ratios 



Tables III gives ratios of branching fractions that show two general patterns. First, the ratios 
are of order unity for modes related by replacing D mesons. The unit ratio for D <^D* replace- 
ment suggests that the additional D* polarization degrees of freedom does not increase produc- 
tion. The unit ratio for •i-^L'W suggests that the production can be described by simple 
isospin relations. Second, the ratios for modes related by addition of pions imply the hierarchy of 

^3-body < •Ss-body < 'S4_body 



V. DECAY DYNAMICS 

A. Three-body decays B^D^*^pp 

Figure 8 gives the Dalitz plot for 3-body B decay modes, where we sum over decays. Each 
point in the vn? {pp)-m? {D^'^ p) plane represents an event and we observe threshold enhancements 



TABLE III: Ratios of branching fractions, r. The (Ttot is the quadrature sum of cXstat and (Tgyst- 
(a) Modes related by s- and d-wave mesons (b) Modes related by charged and neutral D*-*' mesons 



Quantity 



ppTTTT 



/BljOpp 

PpiT 

I ^D^pf-K 

I Bo'^pp-KlX 

IB 



Bd'O 

Bo'+f 

Average with of 8.9/4 



D+f 



0.95±0.13 (0.08) 
1.37±0.17 (0.06) 
1.00±0.13 (0.05) 
0.64±0.19 (0.13) 
1.12±0.25 (0.13) 
1.01±0.07 



ricTtot (cTstat) Quantity 



Bo^ppn /Bd + pptt 

Bd'°pptt /Bd'+pptt 

BD°pp7rir /Bd + 

B £)*Opp^^ I B D*+pp^^ 

Average with of 6.6/3 



r±(Ttot (cTstat) 



1.12±0.12 (0.05) 
0.82±0.12 (0.05) 
1.80±0.44 (0.19) 
1.03±0.28 (0.21) 
1.00±0.08 



(c) Comparison of 3- to 4-body B decay 



(d) Comparison of 4- to 5-body B decay 



Quantity 



Bd*°p 



B 



D^ppiz 



/Bjjopj 



Average with x of 0.5/1 



rzbo-tot (o-stat) Quantity 



3.66±0.48 (0.22) 
3.08±0.62 (0.31) 
3.44±0.38 



PPTT 

/Bd'+pptt 
I Bijapp-^ 

/Bd'Opptt 



Bd+pptttt 

BjJ^ + pp-j^Tl 
B D^ppTTTT 

Bo*^^pp7nT 

Average with of 6.5/3 



r±(Ttot (cTstat) 



0.50±0.10 (0.04) 
0.41±0.06 (0.04) 
0.80±0.14 (0.06) 
0.51±0.14 (0.10) 
0.49±0.05 
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in both variables.^ Red lines are drawn as visual guides: the threshold effect in pp {D^*^^p) is on 
the left (bottom) of the vertical (horizontal) line at 5 (9) GeV^/c^. The two regions do not overlap. 

We investigate the enhancements with mass projection plots using the background subtraction 
technique [43] by weighting pdf components and associated yields, which were described in Eqn. 3 
and given in Table Ila. The weight for the i*^ event, where yi = {rriESi, ^Ei), is 



CssPs{yi) + CsBPB{y^ 
nsPsiVi) +nBPB{yi) 



and (Csa) 



N events 

= E ■ 



Ps{yj)PK{yj 



[nsPsiyj] + nBPB{yj)Y 



(7) 
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(c) B'' D*'^pp in TTiES- A_E signal box 
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(d) — > D*'^pp scaled tties- A_E sideband 



FIG. 8: Dalitz plots oim?{pp) vs. m?{D^p) for two 3-body B decays. Plot in (a) shows 

pp in the m-^a^-A.E signal box and (b) the sideband region with the number of events 
corresponding to the background yield in (a). Corresponding plots for ^ D*^pp are given in 
(c) and (d), respectively. Events live in the shaded allowed kinematic region except for one 
outlier in (d), which failed the re- fit procedure described in footnote d. Red lines are drawn at 
(x,y) of (5 GeV2/c^ 9 GeV2/c^) for (ab) and (5 GeV7c^ 11.5 GeVVc^) for (cd). 



' B candidates are refit with the B-mass constraint so that events lie in the region allowed by kinematic limits. 
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where A = S, B and Csb quantifies the correlation between the signal and background yields. 

Figure 9 shows the differential branching fractions for pp and D^'^p in different regions of the 
complementary variable, which are found by substituting Wi for ns — Uf, in Eqn. 5 and computing 
e for the appropriate region binned in the two variables. In general, the distributions for D^pp 
and D*^pp are similar. One exception is the comparison of D^*^p in the region w?{pp) > 5GeV^/c^ 
shown in Fig. 9b and f, where two broad peaks around lOOMeV/c^ wide are seen for D^pp at 2.9 
and 3.2 GeV/c^, but only one is seen in D*"pp near threshold. 

Lastly, we consider the possibility that the D'^p threshold effect in — > D'^pp may be due to 
one or both of two Ac states observed [44] at 2.94 and 2.88 GeV/c^ with a full width of 18 and 6 
MeV/c^, respectively. These states would span 1-2 bins in Fig. 9b and are not sufficient to describe 
the 150-200 MeV/c^ wide enhancement. 

B. Four-body decays B — yD^ppiv 

Figure 10 shows the differential branching fractions for four mass variables for 4-body B decays: 
pp, D^p, D^p, and p'n~ . In pp, Figs. lOaeim, we see an enhancement production near threshold 
compared to the phase space expectations. In D^p, Figs. lObfjn, we see no obvious hints of a 
narrow exotic 5-quark state around 3.1GeV/c^ [45]. In D^p, Figs. lOegko, the broad threshold 
structure seen 3-body decays seems to be absent, except perhaps in Fig. 10k for B~ ^ D^ppiT~ . 
In pvr~. Figs. lOdhlp, we observe a narrow structure near 1.5GeV/c^, especially in Fig. lOd for 
B^^D^ppTr~, which refer to as X in subsequent discussions. 

To investigate X further, we show opposite-sign and same-sign pir in finer bins of lOMeV/c^ 
in Figs. 11 and 12 for the two neutral and charged B decay modes, respectively. (The detector 
resolution around 1.5GeV/c^ is less than 4MeV/c^.) In Figs, llac and 12ac showing opposite-sign 
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FIG. 9: 3-body B decay mass projection plots for (abed) B^ —>^ D^pp and (efgh) B^ 
(abef) gives plots of m{D^p) in the regions to the left and right of the vertical red line in 
Fig. 8ab; Fig. (cdgh) plots of m{pp) above and below the horizontal red line. The curve represents 
decays following the uniform phase space model. All requirements are given in GeV^/c"^. 
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pvr", the points with error bars show the total number of events in the m-Es-^E signal box, which 
includes of three components: the B signal without X, the B signal with X, and the non-signal 
events. Although we do not have any way a priori to determine the first component {B signal 
without X) a reasonable procedure, which we adopt, is to take the smoothed histogram pdf from 
the same-sign mass combination of pvr" in the same B decay mode written as Pss- This procedure 
has a systematic uncertainty that is difficult to quantify, because the formation of p and p is not 
neccessarily symmetric with respect to the vr" in the B^D^*>ppTr~ decays. The distribution of pvr" 
events as well as the histogram pdf are shown in Figs, llbd and 12bd. The second component (the 
B signal with X) is described by a Breit-Wigner line shape written as P^u)- The third component 
(the non-signal background) can be determined from the m-Es-^E sidebands scaled to the number 
of non-signal events in the signal box. The scaled sideband event samples is shown in Figs. 11-13 
as a grey-shaded histogram and we note that these are smooth across region around 1.5GeV/c^. 

Figures 11a and 11c show the fits of pvr" ■ for 5° ^ D+ppTT- and B° D*+ ppTT , respectively, 
using Phw and Pgs- Pss is relatively flat across X in Figs, lib and provides a good description of 
non-X in 11a. However, Pss from Fig. lid is not as smooth around the signal region and shows 
a dip just below 1.5GeV/c^. We have not yet tried to quantify the uncertainty on this shape. 
The fit results for D~^ppTT~ and D*'^ppTT~ converged to a mean (width) of 1494.4 it 4.1 (51 it 18) 
MeV/c^ and 1500.8 ± 4.4 (43 ± 17) MeV/c^, respectively, where the errors are statistical. We note 
a small excess of events above 1.6GeV/c^ with respect to Pss, but have not attempted to fit it. 
The statistical significance of X in D~^pp7r~ and D*~^ppTr~ is found by comparing the likelihood 
values of our nominal fit to the background-only hypothesis, y^2 {InCh^^ss — ln>Css), which is 8.6 
and 6.9 standard deviations, respectively, under these background assumptions. 

Figures 12a and 12c show the fit to pvr" for B~ ^ D^pp-K" and B~ D*^pp'K~ , respectively. 
The above analysis is repeated, but the signals for X in these modes are not as clear compared 
to the neutral B decays. Moreover, the fit did not converge for B~ — > D^ppir" following the 
procedure, so we re-fit with the X width fixed to the value found for B^ ^ D~^ppTT~ . The fit results 
are consistent with the values found for neutral B decays, but we leave them out in the average. 

Figure 13 shows alternate fits to piT~ in B^^D+ ppTT . The first. Fig. 13a, uses known excited 
nucleon resonances at 1440, 1520, 1535, and 1650 MeV/c^. The Breit-Wigner parameters are 
fixed to their nominal values and are written as Pn*- The fit with Pss and Pat* does not well 
describe the data around 1.5GeV/c^. The second. Fig. 13b, uses Pbw, Pss, and P^'^, where the 
last component is the same-sign p7r~ distribution from pp-K . This is used to determine 

the systematic uncertainty for our choice of the background pdf, which shifted the mean (width) 
by 0.8 (4) MeV/c^. Lastly, for the mass uncertainty, we add another 0.5MeV/c^ to the absolute 
uncertainty due to the variation in the assumed magnetic field and the detector material based on 
the study of [46]. 

Table IV summarizes the results for X, which can be described by a Breit-Wigner with 

mx = (1497.4 ± 3.0 ± 0.9) MeV/c^ (8) 
Fx = { 47 ±12 ±4 )MeV/c^ 

where the errors are statistical and systematic, respectively. The X yields can be compared with 
the total B yields and their ratios range from 5-15%. 

C. Five-body decays B — ^D^'^ppizTZ 

Figure 14 shows the differential branching fractions for a selection of four mass variables for 5-body 
B decays: pp, D^*^p, D^p, and pvr". In contrast to 3- and 4-body distributions, we see do not see 
wide disagreement with phase space expectations. 
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VI. CONCLUSIONS 



Using BABAR's data set of 455 x 10^ BB pairs, we present the observation and study often baryonic 
i?-meson decays to a D^*\ a proton-antiproton pair, and a system of up to two pious, of which 
six {B~ D^ppTT^, B^ L»*°ppvr", B° D^pp-K'n"^, B^ D*^pptt~tt+, B~ D+ppTT-n-, 
B~ D*^ppTT~Ti~) are first observations. The branching fractions for 3- and 5-body decays are 
suppressed compared to 4-body decays with the hierarchy ^3_body < ^Bs-body < 'S4_body Branching 
fraction ratios for modes related by exchange of D mesons are of order unity. 

The kinematic distributions show a number of pecuUar features in the B sample. For 3-body 
decays, non-overlapping threshold enhancements are seen in m{pp) and m{D^'^p). For 4-body 
decays, m{p7r~) distribution shows a narrow peak with mass of (1497.4 it 3.0 it 0.9) MeV/c^ and 
fuh width of (47 ± 12 ± 4) MeV/c^, where the first (second) errors are statistical (systematic). For 
5-body decays, in contrast to 3- and 4-body decays, mass projections are similar with phase space 
expectations. 

All results are preliminary. 
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FIG. 12: Fits of m{piT) in two charged B decays: (ab) B~ ^ D^ppiT~ and (cd) ^ D*^pp'K~ . 
The plotted sample are events in the signal box of rriEs-^E within 2.5a of the mean; the grey 
histograms are the scaled sidebands. The in-set binning is the same as in the main figure. Caveat 
emptor: The fit for (a) has the width fixed to 51MeV/c^ found in Fig. 11a. 
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(a) m{pn^) fit with pdfs from top to bottom: known TV* (b) m{pn^) fit with pdfs from top to bottom: Breit-Wigner 
resonances at 1440, 1520, 1535, 1650 MeV/c^ (dashed, for X (solid, blue); alternate background from m{pn~) in 
blue); background from m{p-K~) in D^pp-K~ (solid, red). D*^pp-K~ (dashed, red); background from m{p-K~) in 

D^ppTT^ (solid, red). 

FIG. 13: Alternate fits of m(p7r~) in W^D+ pp-K . Known N* resonances are used in (a) and an 
alternate background shape from ppn is used in (b) as a systematic study. The in-set 

binning is the same as in the main figure. The grey histograms are the scaled m-Es-^E sidebands. 
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FIG. 14: 5-body B decay differential branching fractions as functions of pp, D^p, D^^p, and pvr" 
in (abed) ^ D^pp7r~7r~^ , (efgh) B^ ^ D*^ppTT~7:~^ , (ijkl) B~ ^ D~^ppTT~iT^ , and (mnop) 
B^ ^ D*^ppTT^'iT~ , respectively. The smooth curve represents decays following the uniform phase 
space model. 
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